Theoretical Analysis Method of Variable Thickness GFRP Tray

Glass-fiber reinforced polymer (GFRP) bars are increasingly widely used in slope support instead of steel bars or steel pipes. GFRP Bars are generally connected with the slope by combining conical nut and tray, but the tray stress still lacks corresponding theoretical calculation and strength verification methods. Therefore, assuming that the tray is an equal thickness thin plate, the internal force distribution of the tray is calculated using the thin plate bending and cavity expansion theory, and compared with the finite element numerical analysis results of the tray. The calculation and analysis show that the elastic theoretical solution of internal force distribution of equal thickness tray is basically the same as the numerical simulation solution of variable thickness tray. The tray loading and free surface are controlled by hoop tensile and radial compressive stress, respectively. The inner wall of the free surface of the tray is the weakest part of the tray, and the ultimate strength of a GFRP tray is 35.81–53.00% of the standard tensile strength of Φ20 GFRP bars by distortion energy density. This theoretical method can be used for stress analysis of variable thickness trays and has played technical support for promoting the application of GFRP bars in slope support.


Introduction
Anchoring technology is widely used in underground engineering support, and materials of anchor solid generally choose reinforcement first. However, in some areas with corrosive solid and groundwater, reinforcement corrosion is increasingly becoming the primary cause affecting the safety and durability of the anchoring system. Due to the behavior of corrosion resistance, suitable cutting, low density, and high tensile strength, glass fiber reinforced polymer (GFRP) bars are first used in the field of mine roadway support instead of reinforcement [1][2][3][4], and then used in slope support [5,6] and floor anti-floating anchorage [7][8][9][10]. To use GFRP bars more reasonably in anchoring technology, domestic and foreign experts and scholars have carried out a large number of experimental and theoretical studies on the material and mechanical properties of GFRP bars, and the results show that GFRP bars can replace reinforcement in slope support anchorage system [11][12][13][14][15][16][17][18][19].
Compared with steel bolts, the GFRP bolt is easy to cut and has no risk. As a support material, it is more environmentally friendly for underground engineering, foundation pit, and slope support, and can solve many disadvantages of traditional steel bolts. However, the top part of a GFRP bolt is often the weakness of anchorage. The tensile strength of a GFRP bar is high, but the shear resistance is poor. The traditional anchorage will produce stress concentration, resulting in the cutting of GFRP bars, and the early failure of anchorage will make the tensile strength of GFRP bars unable to play normally. Therefore, it is very necessary to design reasonable GFRP bar anchorage. After a lot of research and Trays and conical nuts used in this experiment are made of GFRP, and their specifications are shown in Figure 1 and Table 2, respectively. One side of the tray is a plane (after this, referred to as the loading surface), which contacts the solid. The other side is a conical surface (after this, referred to as the free surface), in which six stiffeners with a rib width of 10 mm are uniformly distributed outside. The tray center has a conical round hole (the hole wall later referred to as the inner wall), and the larger aperture is located on the free surface.

GFRP Bar
Φ20 full-threaded GFRP bars produced by a company in Taiyuan, China this experiment, and their parameters are listed in Table 1. Trays and conical nuts used in this experiment are made of GFRP, and their sp are shown in Figure 1 and Table 2, respectively. One side of the tray is a plan referred to as the loading surface), which contacts the solid. The other side is a con (after this, referred to as the free surface), in which six stiffeners with a rib width o uniformly distributed outside. The tray center has a conical round hole (the ho referred to as the inner wall), and the larger aperture is located on the free surface   20 32.378 6 30 The nut center opens a round hole with the thread arranged on the wall. the outer nut wall is hexagonal and the other is conical, and two 3 mm × 30 m symmetrically opened on both sides of the vertebral body. The anchor system top is formed by passing a GFRP bar through the tray center hole and placin tween them. When the tray is stressed, the force is transmitted to the ancho through the nut and then to the soil around the anchor rod body (As shown in  The nut center opens a round hole with the thread arranged on the wall. One end of the outer nut wall is hexagonal and the other is conical, and two 3 mm × 30 mm slits are symmetrically opened on both sides of the vertebral body. The anchor system at the bolt top is formed by passing a GFRP bar through the tray center hole and placing a nut between them. When the tray is stressed, the force is transmitted to the anchor rod body through the nut and then to the soil around the anchor rod body (As shown in Figure 2). the outer nut wall is hexagonal and the other is conical, and two 3 mm × 3 symmetrically opened on both sides of the vertebral body. The anchor sys top is formed by passing a GFRP bar through the tray center hole and pla tween them. When the tray is stressed, the force is transmitted to the an through the nut and then to the soil around the anchor rod body (As show

Tensile Strength Test of GFRP Bars
Three Φ20 GFRP bars with 75 cm long are tested for tensile strength. The ends of each GFRP bar are tied to a steel wire locator and inserted into a 25 cm Φ32 mm × 2 mm steel sleeve tube, then sealed the outer ends of the sleeve tubes. The expansive cement slurry with 4% expansive agent and 0.5 water-cement ratio is injected into the 25 cm high gap between the bar and sleeve tube several times. At the same time, the sleeve tube wall is struck by a iron rod to ensure the slurry is dense and full, as shown in Figure 3. GFRP specimens are cured for 28 days in an environment with a temperature of 20 ± 0.5 • C and humidity of more than 95%. The tensile strength of GFRP bars is tested by WEW-300B hydraulic universal material testing machine (Shanxi Jingchengsheng Engineering Testing Co., Ltd., Thaiyuan, China) with 2 mm/min.

Tensile Strength Test of GFRP Bars
Three Φ20 GFRP bars with 75 cm long are tested for tensile strength. The ends of each GFRP bar are tied to a steel wire locator and inserted into a 25 cm Φ32 mm × 2 mm steel sleeve tube, then sealed the outer ends of the sleeve tubes. The expansive cement slurry with 4% expansive agent and 0.5 water-cement ratio is injected into the 25 cm high gap between the bar and sleeve tube several times. At the same time, the sleeve tube wall is struck by a iron rod to ensure the slurry is dense and full, as shown in Figure 3. GFRP specimens are cured for 28 days in an environment with a temperature of 20 ± 0.5 °C and humidity of more than 95%. The tensile strength of GFRP bars is tested by WEW-300B hydraulic universal material testing machine (Shanxi Jingchengsheng Engineering Testing Co., Ltd., Thaiyuan, China) with 2 mm/min.

Anchoring of GFRP Bars
The ultimate bearing capacity test of GFRP trays is completed in the soil nailing wall support project of the Xiyang Utility Tunnel Project. Pre-drill three holes with 110 mm diameter and 6.5 m deep in slope soil, intercept three holes with a length of 7 m Φ 20 GFRP bars, install a locator every 1 m along the GFRP bars, put Φ20 GFRP bars into the three holes, and expose 0.5 m, respectively. Then, the cement slurry with a water-cement ratio of 0.6 (raw material is PO 42.5 cement) is injected into holes until filled. A square bearing platform shall be made at the orifice, and the surface of the bearing platform shall be perpendicular to the axis of the GFRP bar. The on-site pull-out test shall be carried out after curing 14 days.

Test Device
The loading device is arranged according to Figure 4. Test equipment includes a YCW60B50 through-core jack (Shanxi Jiangong Construction Engineering Testing Co., Ltd., Taiyuan, China) with 20 cm stroke, electric oil pump, pressure gauge, digital dial indicator with 0.01 mm accuracy and 30 mm range, and steel pad with rubber pad. The jack, tray, and nut are assembled and debugged after passing through the GFRP bars anchored on the slope in turn. The test loading is carried out to ensure the good working of

Anchoring of GFRP Bars
The ultimate bearing capacity test of GFRP trays is completed in the soil nailing wall support project of the Xiyang Utility Tunnel Project. Pre-drill three holes with 110 mm diameter and 6.5 m deep in slope soil, intercept three holes with a length of 7 m Φ20 GFRP bars, install a locator every 1 m along the GFRP bars, put Φ20 GFRP bars into the three holes, and expose 0.5 m, respectively. Then, the cement slurry with a water-cement ratio of 0.6 (raw material is PO 42.5 cement) is injected into holes until filled. A square bearing platform shall be made at the orifice, and the surface of the bearing platform shall be perpendicular to the axis of the GFRP bar. The on-site pull-out test shall be carried out after curing 14 days.

Test Device
The loading device is arranged according to Figure 4. Test equipment includes a YCW60B50 through-core jack (Shanxi Jiangong Construction Engineering Testing Co., Ltd., Taiyuan, China) with 20 cm stroke, electric oil pump, pressure gauge, digital dial indicator with 0.01 mm accuracy and 30 mm range, and steel pad with rubber pad. The jack, tray, and nut are assembled and debugged after passing through the GFRP bars anchored on the slope in turn. The test loading is carried out to ensure the good working of loading equipment and instruments and a reliable connection between the trays and GFRP bars. After unloading, start the test.

Test method
The loading process of this test adopted the multi-cycle loading method. The test process is carried out according to anchor tests in a current specification named 'Technical specification for retaining and protection of building foundation excavations (JGJ20-2012)' [35]. Set a load of each level as 6 kN, and record the reading on the dial indicator every 5 min. During the loading process of each level of load, the next level of load can be applied only when the displacement of the anchor head is not greater than 0.1 mm. In case of GFRP bars damage, tray failure, or nut damage during the test, the test shall be stopped.

Results Analysis of GFRP Bars Tensile Strength Test
There is no significant change in GFRP bars at the loading beginning in these tests. As the load increase to about ( is the ultimate tensile load of GFRP bars), the specimen emits a slight and intermittent "hissing" sound, and the surface still has no obvious change. When the load reaches , the specimen makes noise and destroys. The middle part of the specimen gradually uplifts to form a lantern-shaped, and some broken glass fibers are distributed around the uplift in white radial bundles, as shown in Figure 5. The stress-strain curve of GFRP bars consists of three parts: initial gentle section, linear increasing section, and brittle failure section, as shown in Figure 6. The appearance of the initial flat section of the curve may be caused by the slight slip between the specimen

Test Method
The loading process of this test adopted the multi-cycle loading method. The test process is carried out according to anchor tests in a current specification named 'Technical specification for retaining and protection of building foundation excavations (JGJ20-2012)' [35]. Set a load of each level as 6 kN, and record the reading on the dial indicator every 5 min. During the loading process of each level of load, the next level of load can be applied only when the displacement of the anchor head is not greater than 0.1 mm. In case of GFRP bars damage, tray failure, or nut damage during the test, the test shall be stopped.

Results Analysis of GFRP Bars Tensile Strength Test
There is no significant change in GFRP bars at the loading beginning in these tests. As the load increase to about 1 3 F b (F b is the ultimate tensile load of GFRP bars), the specimen emits a slight and intermittent "hissing" sound, and the surface still has no obvious change. When the load reaches F b , the specimen makes noise and destroys. The middle part of the specimen gradually uplifts to form a lantern-shaped, and some broken glass fibers are distributed around the uplift in white radial bundles, as shown in Figure 5.

Test method
The loading process of this test adopted the multi-cycle loading method. The test process is carried out according to anchor tests in a current specification named 'Technical specification for retaining and protection of building foundation excavations (JGJ20-2012)' [35]. Set a load of each level as 6 kN, and record the reading on the dial indicator every 5 min. During the loading process of each level of load, the next level of load can be applied only when the displacement of the anchor head is not greater than 0.1 mm. In case of GFRP bars damage, tray failure, or nut damage during the test, the test shall be stopped.

Results Analysis of GFRP Bars Tensile Strength Test
There is no significant change in GFRP bars at the loading beginning in these tests. As the load increase to about ( is the ultimate tensile load of GFRP bars), the specimen emits a slight and intermittent "hissing" sound, and the surface still has no obvious change. When the load reaches , the specimen makes noise and destroys. The middle part of the specimen gradually uplifts to form a lantern-shaped, and some broken glass fibers are distributed around the uplift in white radial bundles, as shown in Figure 5. The stress-strain curve of GFRP bars consists of three parts: initial gentle section, linear increasing section, and brittle failure section, as shown in Figure 6. The appearance of the initial flat section of the curve may be caused by the slight slip between the specimen chuck and the fixture without preload during the tensile test. As a composite material The stress-strain curve of GFRP bars consists of three parts: initial gentle section, linear increasing section, and brittle failure section, as shown in Figure 6. The appearance of the initial flat section of the curve may be caused by the slight slip between the specimen chuck and the fixture without preload during the tensile test. As a composite material mainly made of brittle glass fiber, GFRP bars directly change from the elastic stage to the brittle failure stage when reaching the load limit. Due to less glass fiber and more matrix resin, the superficial coat of the GFRP bar has a high slip value and cracks first. brittle failure stage when reaching the load limit. Due to less glass fiber and more matrix resin, the superficial coat of the GFRP bar has a high slip value and cracks first. The ultimate tensile strength of three Φ20 GFRP bars is 448.5 MPa, 372.5 MPa, and 353.5 MPa, with an average value of 391.5 MPa. The average elastic modulus of GFRP bars calculated by the linear increase part of the stress-strain curve is 7 GPa, which is quite different from the records of existing literature that will be greater than 40 GPa [6,19]. This difference may be related to the processing technology and glass fiber content of different manufacturers.

Analysis of Nut-Tray Anchorage Test Results
Three groups of the nut-tray anchorage were tested on-site. In the tests, GFRP bars were not damaged, nuts were cracked to varying degrees, and the relative displacements between nuts and trays were about 10 mm. There is extrusion deformation in the inner wall of the loading surface of the tray. The inner wall of the free surface of the tray has obvious cracks, which are distributed radially along the hole wall, as shown in Figure 7. According to the above phenomena, the failure mode of the nut-tray anchorage test is mainly nut cracking and nut sliding with GFRP bars. The damage of the tray is the local cracking of the inner wall, while the tray surface has no obvious damage and serious deformation, and the deformation of the tray surface is a recoverable small deformation. The ultimate tensile strength of three Φ20 GFRP bars is 448.5 MPa, 372.5 MPa, and 353.5 MPa, with an average value of 391.5 MPa. The average elastic modulus of GFRP bars calculated by the linear increase part of the stress-strain curve is 7 GPa, which is quite different from the records of existing literature that will be greater than 40 GPa [6,19]. This difference may be related to the processing technology and glass fiber content of different manufacturers.

Analysis of Nut-Tray Anchorage Test Results
Three groups of the nut-tray anchorage were tested on-site. In the tests, GFRP bars were not damaged, nuts were cracked to varying degrees, and the relative displacements between nuts and trays were about 10 mm. There is extrusion deformation in the inner wall of the loading surface of the tray. The inner wall of the free surface of the tray has obvious cracks, which are distributed radially along the hole wall, as shown in Figure 7. According to the above phenomena, the failure mode of the nut-tray anchorage test is mainly nut cracking and nut sliding with GFRP bars. The damage of the tray is the local cracking of the inner wall, while the tray surface has no obvious damage and serious deformation, and the deformation of the tray surface is a recoverable small deformation.
The tray maximum loading pressures are 90 kN, 96 kN, and 96 kN, with an average value of 94 kN. Therefore, the tensile strength of GFRP bars in this test is between 286.62 MPa and 305.73 MPa, about 73.21-78.09% of the standard tensile strength (391.5 MPa) of Φ20 GFRP bars and does not exceed 85%. After unloading, the relative displacement between the nut and tray is about 10 mm. In the three tests, one tray produces obvious cracks distributed along radial direction near the free surface on the inner wall, and others have no noticeable damage, as shown in Figure 7.

Theoretical Analysis on Tray Stress
In the ultimate bearing capacity test of variable thickness tray, the through center Jack applies a transverse load perpendicular to the loading surface of the tray, and the external force of the tray is transmitted to the GFRP bar through the nut and transformed into the axial load of the GFRP bar. The sidewall of the center hole of the tray forms a circumferential constraint on the nut and bears the extrusion effect caused by the expansion of the nut. In order to theoretically analyze the deformation and internal force of the tray, the variable thickness tray is simplified into an equal thickness circular plate with a central circular hole, and the stress condition of the tray under transverse load is analyzed by thin plate bending theory; The small hole expansion theory is used to analyze the extrusion effect of the nut on the central hole wall of the tray. The internal force and deformation results of the tray can be obtained by superposition of the above two analysis results. Based on the above simplification, it is assumed that the central symmetry axis of the circular plate with equal thickness is the z-axis, and the axis perpendicular to the zaxis along the radial direction of the thickness bisector of the circular plate is defined as ρ-axis, thereby establishing a coordinate system as shown in Figure 8. Through the above analysis, the mechanical analysis of the tray is carried out in the following three steps: first, when calculating the internal force of the tray by using the thin plate bending theory, it is assumed that the tray is integrated with the nut and GFRP bar, and the relative sliding between GFRP bar, nut and tray is not considered; The second step is to simplify the inner wall of the central hole of the tray into an inclined plane, and the included angle between the inclined plane and the z-axis is θ. There are a friction force f parallel to the hole wall and a positive pressure N perpendicular to the hole wall between the nut and the tray hole wall. The friction force f and positive pressure N are decomposed into loads f1 and N1 parallel to the z-axis, respectively, along the radial direction of the tray ρ-axial loads f2

Theoretical Analysis on Tray Stress
In the ultimate bearing capacity test of variable thickness tray, the through center Jack applies a transverse load perpendicular to the loading surface of the tray, and the external force of the tray is transmitted to the GFRP bar through the nut and transformed into the axial load of the GFRP bar. The sidewall of the center hole of the tray forms a circumferential constraint on the nut and bears the extrusion effect caused by the expansion of the nut. In order to theoretically analyze the deformation and internal force of the tray, the variable thickness tray is simplified into an equal thickness circular plate with a central circular hole, and the stress condition of the tray under transverse load is analyzed by thin plate bending theory; The small hole expansion theory is used to analyze the extrusion effect of the nut on the central hole wall of the tray. The internal force and deformation results of the tray can be obtained by superposition of the above two analysis results. Based on the above simplification, it is assumed that the central symmetry axis of the circular plate with equal thickness is the z-axis, and the axis perpendicular to the z-axis along the radial direction of the thickness bisector of the circular plate is defined as ρ-axis, thereby establishing a coordinate system as shown in Figure 8. Through the above analysis, the mechanical analysis of the tray is carried out in the following three steps: first, when calculating the internal force of the tray by using the thin plate bending theory, it is assumed that the tray is integrated with the nut and GFRP bar, and the relative sliding between GFRP bar, nut and tray is not considered; The second step is to simplify the inner wall of the central hole of the tray into an inclined plane, and the included angle between the inclined plane and the z-axis is θ. There are a friction force f parallel to the hole wall and a positive pressure N perpendicular to the hole wall between the nut and the tray hole wall.
The friction force f and positive pressure N are decomposed into loads f 1 and N 1 parallel to the z-axis, respectively, along the radial direction of the tray ρ-axial loads f 2 and N 2 . F sρ is the equivalent aperture shear force, which is used to balance the transverse uniformly distributed load acting on the tray. Load q 2 has an extrusion effect on the hole wall of the tray, and the cavity expansion is used to calculate the internal force of the equal thickness plate; In the third part, the internal force and deformation of the tray are calculated by the above two parts are superimposed to obtain the internal force distribution and deformation of the tray. Thus: where ϕ is the internal friction angle between the nut and the inner wall, tan ϕ = 0.3. θ is inner wall inclination angle with a value of 0.06922.
wall of the tray, and the cavity expansion is used to calculate the internal force of the equal thickness plate; In the third part, the internal force and deformation of the tray are calculated by the above two parts are superimposed to obtain the internal force distribution and deformation of the tray. Thus: where φ is the internal friction angle between the nut and the inner wall, tanφ = 0.3. θ is inner wall inclination angle with a value of 0.06922.

Force Analysis of Tray Bending under Transverse Load
The transverse load exerted by the jack on the tray loading surface is simplified as a uniform load, and the contact between the nut and the tray inner wall is facilitated as a fixed bearing. According to the field test results, the tray ultimate failure load is 94 kN, the elastic modulus is 7 GPa, and the Poisson's ratio is 0.2. The tray can be simplified into two types of equal thickness thin plates. The thickness of thin plates A and B select the tray edge height as 12.16 mm and the tray inner wall height as 35.34 mm, respectively. The inner diameter of the tray is 14.4 mm, and the outer is 69.9 mm. According to the tray loading surface area, the equivalent uniform load on the tray loading surface is 6.40 MPa. According to reference [36], the bending deflection of an axisymmetric circular thin plate can represent Equation (4): where C1, C2, C3, and C4 are undetermined coefficients. ω is the bending deflection, m. ρ is the horizontal distance between the tray and the symmetric axis, m. D is bending stiffness, kN·m. From Equation (4), the equations of rotation, bending moment, and shear force of the tray along the radial direction are = C + 2ρlnρC + ρC + 2ρC + ρ (5)

Force Analysis of Tray Bending under Transverse Load
The transverse load exerted by the jack on the tray loading surface is simplified as a uniform load, and the contact between the nut and the tray inner wall is facilitated as a fixed bearing. According to the field test results, the tray ultimate failure load is 94 kN, the elastic modulus is 7 GPa, and the Poisson's ratio is 0.2. The tray can be simplified into two types of equal thickness thin plates. The thickness of thin plates A and B select the tray edge height as 12.16 mm and the tray inner wall height as 35.34 mm, respectively. The inner diameter of the tray is 14.4 mm, and the outer is 69.9 mm. According to the tray loading surface area, the equivalent uniform load q 1 on the tray loading surface is 6.40 MPa.
According to reference [36], the bending deflection of an axisymmetric circular thin plate can represent Equation (4): where C 1 , C 2 , C 3 , and C 4 are undetermined coefficients. ω is the bending deflection, m. ρ is the horizontal distance between the tray and the symmetric axis, m. D is bending stiffness, kN·m. From Equation (4), the equations of rotation, bending moment, and shear force of the tray along the radial direction are where M ρ is the radial bending moment of the tray, kN·m/m. M ϕ is the hoop bending moment of the tray, kN·m/m. F sρ is the tray radial shear force, kN/m. µ is Poisson's ratio. E is elastic modulus, MPa. δ is thin plate thickness, m. The boundary conditions of the tray are: According to the above equations and boundary conditions, the deflection equations of thin plates A and B can be obtained, respectively: The radial deflection curves of the tray are shown in Figure 9. Under the same load conditions and boundary conditions, the tray thickness has a significant influence on the radial deflection amplitude of the tray.
where M is the radial bending moment of the tray, kN·m/m. M is the hoop bendin moment of the tray, kN·m/m. F is the tray radial shear force, kN/m. μ is Poisson's ratio E is elastic modulus, MPa. δ is thin plate thickness, m. The boundary conditions of the tray are: According to the above equations and boundary conditions, the deflection equation of thin plates A and B can be obtained, respectively: ω = −2.533 × 10 lnρ − 3.575ρ lnρ − 7.301ρ − 1.237 × 10 + 91.460ρ (12 The radial deflection curves of the tray are shown in Figure 9. Under the same load conditions and boundary conditions, the tray thickness has a significant influence on th radial deflection amplitude of the tray. The radial stress and hoop stress of the tray are: When the thickness is reduced by 2.9 times from 35.34 mm to 12.16 mm, the deflection at the tray edge increases by 24.5 times from 0.30 mm to 7.35 mm.
The radial stress and hoop stress of the tray are: where σ ρ is the radial stress of the tray, MPa. σ ϕ is the hoop stress of the tray, MPa. Referring to Figure 10, the radial stress of the tray loading surface is tensile and gradually decreases from the center to the tray edge. At 58 mm away from the symmetrical axis, the radial stresses of thin plates A and B are close to 0. The radial stress attenuation amplitude of thin plate A is faster than B, and the maximum radial stress at the inner wall of thin plate A and B is 773.04 MPa and 91.52 MPa, respectively, with a difference of nearly 8.5 times.
where σ is the radial stress of the tray, MPa. σ is the hoop stress of the tray, MPa. Referring to Figure 10, the radial stress of the tray loading surface is tensile and gradually decreases from the center to the tray edge. At 58 mm away from the symmetrical axis, th radial stresses of thin plates A and B are close to 0. The radial stress attenuation amplitude o thin plate A is faster than B, and the maximum radial stress at the inner wall of thin plate A and B is 773.04 MPa and 91.52 MPa, respectively, with a difference of nearly 8.5 times. The hoop stress of the free surface of thin plates A and B is tensile stress, and th variation law is the same. The tensile from the center to the edge of the thin plate firs increases and then decreases. The difference is that the hoop stress of thin plate A i greater than that of thin plate B. The hoop stresses of thin plate A and B reach the peak a 23 mm from the symmetry axis, which are 248. 30  Calculate the radial shear F of the pallet from Equation (8), The results are plotted in Figure 11. The tray thickness does not affect the tray radial shear force. The maximum shea force is 1038.93 kN/m at the inner wall and gradually decreases to 0 along the radial direction The hoop stress of the free surface of thin plates A and B is tensile stress, and the variation law is the same. The tensile from the center to the edge of the thin plate first increases and then decreases. The difference is that the hoop stress of thin plate A is greater than that of thin plate B. The hoop stresses of thin plate A and B reach the peak at 23 mm from the symmetry axis, which are 248. 30  Calculate the radial shear F sρ of the pallet from Equation (8), The results are plotted in Figure 11. The tray thickness does not affect the tray radial shear force. The maximum shear force is 1038.93 kN/m at the inner wall and gradually decreases to 0 along the radial direction.
The resultant shear stress of the tray inner wall Calculate the radial shear F of the pallet from Equation (8), The results are plotted i Figure 11. The tray thickness does not affect the tray radial shear force. The maximum shea force is 1038.93 kN/m at the inner wall and gradually decreases to 0 along the radial direction

Influence of Nut Extrusion on the Tray Internal Force
The extrusion effect of the nut on the inner wall of the tray can be simplified as the cavity expansion problem. From the previous analysis, the squeezing effect on the inner wall of the tray comes from two aspects: the radial component generated by the positive pressure perpendicular to the inner wall of the tray and the radial component caused by the friction parallel to the inner wall of the tray. The resultant force of these two components jointly produces radial expansion on the tray, as shown in Figure 8.
There is friction contact between the nut and the tray inner wall. Both the nut and the tray are GFRP composites, whose friction coefficient is mainly related to the matrix material, glass fiber content, and stress conditions. According to the existing literature [37][38][39][40], although affected by glass fiber content and test conditions, the friction coefficient has been between 0.27 and 0.34. Therefore, in the analysis, the friction coefficient between the nut and the tray can be selected as 0. 30. In the calculation, the tray inner radius is determined according to the equivalent aperture, which is selected as the average value ρ 1 = 15.625 mm. The hole height (h) is determined according to the actual height of the tray inner wall, h = 35.34 mm. According to the unfavorable factors, the influence range of small hole expansion is the radius of the small circular surface of the tray, ρ 2 = 32.33 mm. Therefore, from Equation (1) to Equation (3), the radial stress q 2 acted on the equivalent radius inner wall can be calculated as 71.83 MPa.
The tray stress under extrusion is calculated according to Equations (16) and (17), and its distribution is shown in Figure 12: where ρ 1 is the tray equivalent inner diameter, m. ρ 2 is the maximum equivalent radius of the cavity expansion effect, m. ρ is the distance between the calculated point and the tray symmetrical axis, m. σ ϕ is the hoop stress produced by the cavity expansion, MPa. σ ρ is the radial stress created by the cavity expansion, MPa. q 2 is the radial load applied to the tray inner wall, MPa. The radial and hoop stress caused by the nut extrusion on the tray is compressive and tensile, respectively, as shown in Figure 12. The maximum hoop and radial stresses are 115.62 MPa and 71.83 MPa, respectively, which are at the tray inner wall, and then gradually decrease along the tray radial direction. The hoop stress plays a leading role along the tray radial direction and is larger than the radial stress. At the equivalent radius of 32.33 mm, the radial stress is 0, and the hoop tensile is 43.78 MPa.

Superposition of Tray Internal Force
From the previous analysis, the internal force distribution of the tray reflects the extrusion effect of the nut on the inner wall of the tray and the bending effect of the tray surface under the transverse uniformly distributed load. Therefore, the internal forces generated by them should be superimposed during radial and hoop stresses of the tray are calculated by theory. According to the cavity expansion results, the stresses produced by nut extrusion on the tray are radial compressive and hoop tensile stress. For transverse load, the radial and hoop stresses of the tray loading and free surface are tensile and compressive stress, respectively. Stresses on two surfaces are mirror symmetry according to the transverse axis. The results of superposition of thin plate bending and small hole expansion are as follows: for the stress surface, the hoop tensile stress increases and the radial stress decreases; For the free surface, on the contrary, the hoop stress decreases, and the radial stress increases.
The tray stress superposition analysis of radial and hoop stresses is carried out with the thin plate B as an example, as shown in Figure 13a. For the radial stress of the tray loading surface, on the one hand, the radial compressive stress caused by the cavity expansion becomes 0 at the distance of 32.33 mm from the tray symmetry axis, and its influence range is small and numerical attenuation is fast. After superimposing the radial tensile stress caused by thin plate bending and the radial compressive stress caused by cavity expansion, the resultant radial stress is tensile. On the other hand, although the radial compressive stress caused by cavity expansion has a minor influence range, its stress value near the tray inner wall is almost close to the radial tensile stress caused by the thin plate bending. Therefore, the superposition makes the radial stress smaller, which is beneficial to the tray inner wall. The trend of resultant radial stress increases first and then decreases. The pressure at the tray inner wall has a significant influence, which is 8.51 MPa. The resultant maximum radial stress is only 20.27 MPa. The radial and hoop stress caused by the nut extrusion on the tray is compressive and tensile, respectively, as shown in Figure 12. The maximum hoop and radial stresses are 115.62 MPa and 71.83 MPa, respectively, which are at the tray inner wall, and then gradually decrease along the tray radial direction. The hoop stress plays a leading role along the tray radial direction and is larger than the radial stress. At the equivalent radius of 32.33 mm, the radial stress is 0, and the hoop tensile is 43.78 MPa.

Superposition of Tray Internal Force
From the previous analysis, the internal force distribution of the tray reflects the extrusion effect of the nut on the inner wall of the tray and the bending effect of the tray surface under the transverse uniformly distributed load. Therefore, the internal forces generated by them should be superimposed during radial and hoop stresses of the tray are calculated by theory. According to the cavity expansion results, the stresses produced by nut extrusion on the tray are radial compressive and hoop tensile stress. For transverse load, the radial and hoop stresses of the tray loading and free surface are tensile and compressive stress, respectively. Stresses on two surfaces are mirror symmetry according to the transverse axis. The results of superposition of thin plate bending and small hole expansion are as follows: for the stress surface, the hoop tensile stress increases and the radial stress decreases; For the free surface, on the contrary, the hoop stress decreases, and the radial stress increases.
The tray stress superposition analysis of radial and hoop stresses is carried out with the thin plate B as an example, as shown in Figure 13a. For the radial stress of the tray loading surface, on the one hand, the radial compressive stress caused by the cavity expansion becomes 0 at the distance of 32.33 mm from the tray symmetry axis, and its influence range is small and numerical attenuation is fast. After superimposing the radial tensile stress caused by thin plate bending and the radial compressive stress caused by cavity expansion, the resultant radial stress is tensile. On the other hand, although the radial compressive stress caused by cavity expansion has a minor influence range, its stress value near the tray inner wall is almost close to the radial tensile stress caused by the thin plate bending. Therefore, the superposition makes the radial stress smaller, which is beneficial to the tray inner wall. The trend of resultant radial stress increases first and then decreases. The pressure at the tray inner wall has a significant influence, which is 8.  For the hoop stress of the loading surface of the tray, the stresses generated by the thin plate bending and the cavity expansion are all tensile. The resultant hoop stress increase significantly, its trend decreases from the tray inner wall to the outside, and the tray equivalent inner wall (15.625 mm from the symmetrical axis) has the largest stress (138.09 MPa). Within the range of stress superposition, the contribution rate of transverse load in the resultant hoop tensile stress after superposition is 16~37%, and the contribution rate of cavity expansion is 8~63%. Among them, the difference in the contribution rate at the inner wall of the tray is the largest. The contribution rates of cavity expansion and transverse load are 84% and 16%, respectively. Cavity expansion plays a leading role in the hoop stress of the tray and cannot be ignored.
The same method is used to analyze the free surface combined stress of the thin plate B. The radial stresses produced by the thin plate bending and the cavity expansion are all compressive stress, and the resultant radial stress decreases from the tray inner wall to the outside. The tray inner wall (16.85 mm from the symmetrical axis) has the maximum radial stress of 128.37 MPa. The proportion of the stresses produced by the thin plate bending and the cavity expansion is about 55% and 45%, respectively. The properties of hoop stress caused by transverse load and cavity expansion are just opposite, one is compressive stress and the other is tensile stress. The compressive stress generated by thin plate bending has a weakening effect on the tensile stress caused by cavity expansion, which is reduced by about 33 %. The resultant hoop stress of the free surface is tensile stress, which decreases from the tray inner wall to the outside. The maximum stress appears in the tray inner wall (16.85 mm from symmetrical axis, 75.77 MPa). The radial compressive stress of the tray free surface is much larger than the hoop tensile stress, about 1.69 times, so the radial stress controls the free surface, as shown in Figure 13b.
In summary, the cavity expansion effect is significant on the loading surface and free surface of the tray, so it cannot be ignored. The hoop tensile stress of the loading surface is greater than the radial tensile stress, and the radial compressive stress of the free surface is greater than the hoop tensile stress. Therefore, the loading surface of the thin plate B is controlled by the hoop tensile stress. On the contrary, the free surface of the thin plate B is controlled by the radial compressive stress controls. Comparison from the numerical absolute value, the hoop stress on the loading surface of thin plate B is higher than the radial stress on the free surface. We can also draw a similar conclusion for thin plate A.

Tray Force Analysis of the Numerical Simulation
GFRP tray has a thick middle and thin edge, which is not an ideal thin plate strictly. To analyze the tray limit state stress, Ansys Workbench 2020 R2 software is used to simulate the deformation and internal force of the tray under ultimate transverse load. Firstly, For the hoop stress of the loading surface of the tray, the stresses generated by the thin plate bending and the cavity expansion are all tensile. The resultant hoop stress increase significantly, its trend decreases from the tray inner wall to the outside, and the tray equivalent inner wall (15.625 mm from the symmetrical axis) has the largest stress (138.09 MPa). Within the range of stress superposition, the contribution rate of transverse load in the resultant hoop tensile stress after superposition is 16~37%, and the contribution rate of cavity expansion is 8~63%. Among them, the difference in the contribution rate at the inner wall of the tray is the largest. The contribution rates of cavity expansion and transverse load are 84% and 16%, respectively. Cavity expansion plays a leading role in the hoop stress of the tray and cannot be ignored.
The same method is used to analyze the free surface combined stress of the thin plate B. The radial stresses produced by the thin plate bending and the cavity expansion are all compressive stress, and the resultant radial stress decreases from the tray inner wall to the outside. The tray inner wall (16.85 mm from the symmetrical axis) has the maximum radial stress of 128.37 MPa. The proportion of the stresses produced by the thin plate bending and the cavity expansion is about 55% and 45%, respectively. The properties of hoop stress caused by transverse load and cavity expansion are just opposite, one is compressive stress and the other is tensile stress. The compressive stress generated by thin plate bending has a weakening effect on the tensile stress caused by cavity expansion, which is reduced by about 33 %. The resultant hoop stress of the free surface is tensile stress, which decreases from the tray inner wall to the outside. The maximum stress appears in the tray inner wall (16.85 mm from symmetrical axis, 75.77 MPa). The radial compressive stress of the tray free surface is much larger than the hoop tensile stress, about 1.69 times, so the radial stress controls the free surface, as shown in Figure 13b.
In summary, the cavity expansion effect is significant on the loading surface and free surface of the tray, so it cannot be ignored. The hoop tensile stress of the loading surface is greater than the radial tensile stress, and the radial compressive stress of the free surface is greater than the hoop tensile stress. Therefore, the loading surface of the thin plate B is controlled by the hoop tensile stress. On the contrary, the free surface of the thin plate B is controlled by the radial compressive stress controls. Comparison from the numerical absolute value, the hoop stress on the loading surface of thin plate B is higher than the radial stress on the free surface. We can also draw a similar conclusion for thin plate A.

Tray Force Analysis of the Numerical Simulation
GFRP tray has a thick middle and thin edge, which is not an ideal thin plate strictly. To analyze the tray limit state stress, Ansys Workbench 2020 R2 software is used to simulate the deformation and internal force of the tray under ultimate transverse load. Firstly, the theoretical calculation and numerical simulation of thin plates A and B under transverse load are compared to verify the rationality of the tray parameter selection. Then according to the geometric size of Figure 1, the axisymmetric numerical model is established to simulate the tray.
The previous Section 3. 2 show that under the maximum load of the test, the inner wall of the tray appears cracks, and the tray surface is not obviously damaged. The deformation of the tray surface is recoverable. It can be considered that the tray is still in the elastic stage under the maximum load. Therefore, the linear static analysis is used in the numerical simulation of the tray force in this section. In the numerical simulation, the elastic modulus was measured by Section 3.1, that is, the elastic modulus is 7.0 GPa. The parameters of the tray, nut, and GFRP bar, such as the geometric size, the spatial contact relationship, the boundary conditions, and the applied load, are the same as those of Section 3.2.

Rationality Verification of the Numerical Simulation Parameters
The calculation results of the loading surface of thin plates A and B are shown in Figure 14. The numerical simulation of radial and hoop stress of thin plates A and B are consistent with the theoretical calculation, which is different 0.48~1.77%. Therefore, the numerical model element set and the material parameter value are appropriate. the theoretical calculation and numerical simulation of thin plates A and B under transverse load are compared to verify the rationality of the tray parameter selection. Then according to the geometric size of Figure 1, the axisymmetric numerical model is established to simulate the tray. The previous Section 3. 2 show that under the maximum load of the test, the inner wall of the tray appears cracks, and the tray surface is not obviously damaged. The deformation of the tray surface is recoverable. It can be considered that the tray is still in the elastic stage under the maximum load. Therefore, the linear static analysis is used in the numerical simulation of the tray force in this section. In the numerical simulation, the elastic modulus was measured by Section 3.1, that is, the elastic modulus is 7.0 GPa. The parameters of the tray, nut, and GFRP bar, such as the geometric size, the spatial contact relationship, the boundary conditions, and the applied load, are the same as those of Section 3.2.

Rationality Verification of the Numerical Simulation Parameters
The calculation results of the loading surface of thin plates A and B are shown in Figure 14. The numerical simulation of radial and hoop stress of thin plates A and B are consistent with the theoretical calculation, which is different 0.48~1.77%. Therefore, the numerical model element set and the material parameter value are appropriate.

Numerical Simulation Results of the Tray
The numerical model carries out solid 3D modeling according to the actual size (as shown in Figure 1) and contact relationship of the tray. A nut is set between the tray and GFRP bar, and a thread is set between the nut and GFRP reinforcement. The friction contact element is set between the tray and the nut, in which the contact stiffness is 1 and the friction coefficient is 0.3. The element of the numerical model adopts tetrahedral solid element solid187. The numerical model is divided into 84548 grid elements, and the places with large stress and sudden change are locally densified. The density is 2.0 g/mm and Poisson's ratio is 0.2. The stress-strain relationship adopts the test results of the GFRP pullout test, and the elastic modulus is 7 GPa. The boundary condition is that one end of the GFRP bar is fixed displacement and the other end is a free boundary. The loading method adopts one-time loading, and the uniformly distributed load is applied on the loading surface of the tray, with the size of 6.40 MPa. The simulation time is 195.19 s.
The tray radial deflection curve obtained by numerical simulation is shown in the natural line in Figure 9, between the theory deflection of the thin plate A and B. The tray deflection increases gradually from the inner wall to the outside, and the outer edge has

Numerical Simulation Results of the Tray
The numerical model carries out solid 3D modeling according to the actual size (as shown in Figure 1) and contact relationship of the tray. A nut is set between the tray and GFRP bar, and a thread is set between the nut and GFRP reinforcement. The friction contact element is set between the tray and the nut, in which the contact stiffness is 1 and the friction coefficient is 0.3. The element of the numerical model adopts tetrahedral solid element solid187. The numerical model is divided into 84548 grid elements, and the places with large stress and sudden change are locally densified. The density is 2.0 g/mm and Poisson's ratio is 0.2. The stress-strain relationship adopts the test results of the GFRP pull-out test, and the elastic modulus is 7 GPa. The boundary condition is that one end of the GFRP bar is fixed displacement and the other end is a free boundary. The loading method adopts one-time loading, and the uniformly distributed load is applied on the loading surface of the tray, with the size of 6.40 MPa. The simulation time is 195.19 s.
The tray radial deflection curve obtained by numerical simulation is shown in the natural line in Figure 9, between the theory deflection of the thin plate A and B. The tray deflection increases gradually from the inner wall to the outside, and the outer edge has the maximum deflection (1.17 mm). The tray edge simulation value is 0.16 and 3.9 times the theoretical values of the thin plate A (0.30 mm) and B (7.35 mm). Therefore, the deflection simulation value of the tray is closer to the thin plate B along the radial direction. The initial section of the slope of the tray deflection simulation value is small, which increases significantly between 35 mm and 50 mm from the symmetry axis. The variation range of the deflection simulation value slope is basically consistent with the range of the tray thickness thinning area (32.33-44.40 mm). Therefore, the variation of the tray radial thickness has a significant impact on the radial deflection.
The numerical simulation results of the radial stress and hoop stress of the tray are shown in Figure 15a,b, and the radial stress and hoop stress of the corresponding loading surface and free surface are shown in Figure 15c. It can be seen from Figure 15 that the hoop stress on the loading surface is mainly tensile. From the inner wall (14.4 mm away from the symmetric axis) to the edge of the tray, the hoop stress gradually decreases, and the maximum tensile of the inner wall is 126.43 MPa.  1.17 mm). The tray edge simulation value is 0.16 and 3.9 times the theoretical values of the thin plate A (0.30 mm) and B (7.35 mm). Therefore, the deflection simulation value of the tray is closer to the thin plate B along the radial direction. The initial section of the slope of the tray deflection simulation value is small, which increases significantly between 35 mm and 50 mm from the symmetry axis. The variation range of the deflection simulation value slope is basically consistent with the range of the tray thickness thinning area (32.33-44.40 mm). Therefore, the variation of the tray radial thickness has a significant impact on the radial deflection.
The numerical simulation results of the radial stress and hoop stress of the tray are shown in Figure 15a,b, and the radial stress and hoop stress of the corresponding loading surface and free surface are shown in Figure 15c. It can be seen from Figure 15 that the hoop stress on the loading surface is mainly tensile. From the inner wall (14.4 mm away from the symmetric axis) to the edge of the tray, the hoop stress gradually decreases, and the maximum tensile of the inner wall is 126.43 MPa.  The radial stress of the tray loading surface is mainly tensile and becomes compressive stress after about 62 mm from the symmetric axis, which value increases first and then decreases. The maximum tensile stress (54.983 MPa) is near 42 mm from the symmetry axis, and the maximum compressive stress (1.24 MPa) is near the outer tray edge. Radial compressive stress dominates in the free surface and has a bimodal curve. The first stress peak (146.44 MPa) is at the tray inner wall (16.85 mm away from the symmetric axis), and the second (120.57 MPa) is at 44.40 mm away from the symmetric axis. From the tray inner wall to the outside, the hoop stress of the free surface changes from tensile stress to compressive and then to tensile. The maximum tensile stress at the center hole wall of the tray is 23.67 MPa, and the maximum compressive stress at 44.40 mm away from the symmetry axis is 26.81 MPa. The maximum stress of the tray is the radial compressive stress at the junction of the inner wall and the free surface, with the size of 146.44 MPa. The maximum stress is consistent with the location of the tray crack, as shown in Figure 7. Figure 1 shows that 44.40 mm away from the symmetry axis is the thickness turning area of the tray free surface from thick to thin (after this called thickness turning zone, TTZ), and thus the tray center hole and the thickness turning area has an important influence on the tray stress distribution. As known in the stress envelope, respectively, taking the hoop tensile and radial compressive stress distribution line of the loading and free surface as the boundary, the tray stress analysis can control all unfavorable factors.
The radial and hoop stresses of the tray inner wall are shown in Figure 16. The hoop and radial stresses of the tray inner wall are tensile and compressive stress, respectively, which are consistent with the theoretical analysis result of the cavity expansion. The distribution of hoop tensile and radial compressive stress is opposite. From the loading surface of the tray to the free surface, the hoop tensile stress decreases from 126.43 MPa to 23.67 MPa, 85.44 MPa. The radial stress is generally compressive, and the local tensile stress is near the loading surface. The maximum tensile stress is 4.14 MPa, and the maximum compressive stress is 146.44 MPa near the free surface, with an average of 67.62 MPa. Compared with the theoretical calculation of cavity expansion, the hoop tensile stress theoretical value of the inner wall is 115.51 MPa, and the error is 26.03%; The theoretical value of the radial compressive stress of the inner wall is 71.83 MPa, and the error is 5.87%. The theoretical solution of the radial stress on the tray inner wall is close to the numerical simulation results, however, the hoop stress isn't. The reasons are the following two points. (1) The center hole and the radial thickness of the tray are changed in the numerical simulation but use the equivalent radius in the theoretical calculation. Therefore, there is a particular deviation. (2) The influence of the wall thickness on the hoop tensile stress is more considerable than radial compressive stress. When the distance from the inner wall is close to the hole wall thickness, the radial stress tends to 0. however, the hoop stress is still considerable. The wall thickness is changed in the numerical analysis of the tray. Compared with the wall thickness of the theoretical calculation, the wall thickness of the loading surface changes significantly, and that of the free surface changes little. Therefore, the hoop stress deviation is relatively large while the radial is relatively small.

Comparison of Theoretical Calculation and Numerical Simulation
The similarities and differences between theoretical analysis and numerical simulation in terms of geometric parameters, physical parameters, boundary conditions, load type and size, analysis method, and calculation results are shown in Table 3. The elastic modulus, Poisson's ratio, and friction coefficient of the two analysis methods are the same, the boundary conditions and load types and sizes are the same, and the analysis methods are linear elasticity. However, the geometric parameters of the tray used in the two analysis methods are quite different. The geometric model used in the numerical simulation is modeled according to the actual size of the tray. The radial thickness of the tray is variable thickness, and the central opening is conical. In order to adopt elastic mechanics analysis, the thickness of the tray is simplified to equal thickness in theoretical analysis, and the central opening is simplified to an equal diameter circular hole.

Comparison of Theoretical Calculation and Numerical Simulation
The similarities and differences between theoretical analysis and numerical simulation in terms of geometric parameters, physical parameters, boundary conditions, load type and size, analysis method, and calculation results are shown in Table 3. The elastic modulus, Poisson's ratio, and friction coefficient of the two analysis methods are the same, the boundary conditions and load types and sizes are the same, and the analysis methods are linear elasticity. However, the geometric parameters of the tray used in the two analysis methods are quite different. The geometric model used in the numerical simulation is modeled according to the actual size of the tray. The radial thickness of the tray is variable thickness, and the central opening is conical. In order to adopt elastic mechanics analysis, the thickness of the tray is simplified to equal thickness in theoretical analysis, and the central opening is simplified to an equal diameter circular hole.   Comparative conclusions (1) The radial deflection, radial stress, and hoop stress calculated by equal thickness plate B are the same as the numerical simulation results. (2) The tensile stress of the inner wall of the free surface of the tray is the largest, which is the weakest part of the tray, and the above analysis conclusion is the same as the damage location of the tray measured in the nut-tray anchorage test. For the tray free surface, the stress concentration occurs at TTZ in the numerical simulation results, which leads to peak radial compressive stress and the bimodal curve of the compressive stress of the free surface. The theoretical calculation results of the thin plates A and B are all the maximum values at the tray inner wall, which causes the difference between the theoretical calculation and the numerical simulation curves of the radial stress of the free surface. For thin plate A, the radial stress (657.69 MPa) at the tray inner wall (16.85 mm away from the symmetrical axis) is far more than the ultimate tensile strength (391.5 MPa) of GFRP bars measured by the tensile strength test. Therefore, the theoretical calculation ultimate tensile strength of thin plate A is quite different from the test results, that is, the thickness assumption of thin plate A is unreasonable.
For thin plate B, the simulated and theoretical values of radial stress on the free surface at the inner wall of the tray are 146.44 MPa and 128.37 MPa, respectively. The theoretical calculation of thin plate B is 12.34% less than the simulated value. At the tray TTZ, the numerical simulation value (120.57 MPa) of radial stress is pretty different from the theoretical calculation value (5.69 MPa). Therefore, it is suggested that measures should be taken to reduce local stress concentration in the tray design. By comparing the hoop stress of the loading surface and the radial stress of the free surface of the tray, their numerical simulation values are close to the theoretical calculation value of the thin plate B. For the hoop tensile stress of the loading surface, the numerical simulation result is consistent with the theoretical calculation result in the numerical value and the variation law along the radial direction. Without considering the abrupt change of radial compressive stress caused by the drastic shift in tray thickness, at the tray inner wall of thin plate B, the numerical simulation and the theoretical calculation values of radial compressive stress have similar values and distribution laws. Therefore, for the variable thickness tray of Both theoretical calculation and numerical simulation show that the failure is controlled by the most unfavorable factors, in which the loading surface of the tray is controlled by hoop stress and the free surface is controlled by radial stress, as shown in Figure 17. For the tray free surface, the stress concentration occurs at TTZ in the numerical simulation results, which leads to peak radial compressive stress and the bimodal curve of the compressive stress of the free surface. The theoretical calculation results of the thin plates A and B are all the maximum values at the tray inner wall, which causes the difference between the theoretical calculation and the numerical simulation curves of the radial stress of the free surface. For thin plate A, the radial stress (657.69 MPa) at the tray inner wall (16.85 mm away from the symmetrical axis) is far more than the ultimate tensile strength (391.5 MPa) of GFRP bars measured by the tensile strength test. Therefore, the theoretical calculation ultimate tensile strength of thin plate A is quite different from the test results, that is, the thickness assumption of thin plate A is unreasonable.
For thin plate B, the simulated and theoretical values of radial stress on the free surface at the inner wall of the tray are 146.44 MPa and 128.37 MPa, respectively. The theoretical calculation of thin plate B is 12.34% less than the simulated value. At the tray TTZ, the numerical simulation value (120.57 MPa) of radial stress is pretty different from the theoretical calculation value (5.69 MPa). Therefore, it is suggested that measures should be taken to reduce local stress concentration in the tray design. By comparing the hoop stress of the loading surface and the radial stress of the free surface of the tray, their numerical simulation values are close to the theoretical calculation value of the thin plate B. For the hoop tensile stress of the loading surface, the numerical simulation result is consistent with the theoretical calculation result in the numerical value and the variation law along the radial direction. Without considering the abrupt change of radial compressive stress caused by the drastic shift in tray thickness, at the tray inner wall of thin plate B, the numerical simulation and the theoretical calculation values of radial compressive stress have similar values and distribution laws. Therefore, for the variable thickness tray of GFRP, the tray is simplified into a thin plate with equal thickness according to the thickness of the central opening, and the internal force superposition analysis is carried out according to the axisymmetric thin plate bending theory and cavity expansion theory of elasticity, which can theoretically predict the internal force distribution and the main controlling factors of failure of the tray.
It can be seen that the radial deflection, radial stress, and hoop stress of equal thickness plate B calculated by the theoretical analysis have the same variation as the numerical simulation results, and the numerical values are also consistent. The local difference between the theoretical analysis and the numerical simulation results is mainly caused by the difference in the geometric size of the tray.

Ultimate Strength of Tray Failure
According to the tray ultimate bearing capacity test, the maximum loading average value of the tray is 94 kN, which is converted to the GFRP bar strength of 299.22 MPa. The loading of the ultimate bearing capacity test is equivalent to 76.43% of the GFRP bar standard value (391.5 MPa) of the ultimate tensile strength. Considering the discreteness and brittle failure of the tensile strength test value of the GFRP bar, according to the most unfavorable consideration, 84.64% of the minimum test value (353.5 MPa) has been loaded in the tensile strength test, which meets the maximum loading requirement of 85% not exceeding the standard value of the tensile strength specified in the specification JGJ120-2012.
The theoretical calculation and numerical simulation results show that the inner wall and the TTZ of the tray loading and free surface are the places where the stress is enormous. Therefore, the theoretical equation of the distortion energy density, namely Equation (18), can be used to calculate the material strength: where σ 1 , σ 2 and σ 3 are the maximum, intermediate, and minimum principal stresses, respectively, MPa. Calculation results of tray strength are shown in Table 4.  35.81% to 53.00% (average proportion is 44.17%) of 391.5 MPa, these strengths are far less than the Φ20 GFRP bar standard value of ultimate tensile strength. However, the on-site test of the tray shows that there are light cracks on the surface of the tray, which proves that the tray has reached the limit state. Therefore, it is speculated that the ultimate strength of the tray may be inconsistent with that of the GFRP bars material, and the ultimate strength of the tray should be less than the uniaxial tensile strength of the GFRP bars material. The reasons are as follows: (1) The glass fiber distribution in the tray is different from GFRP bars. The glass fiber can be approximately unidirectional in GFRP bars but at least two-dimensional distribution in the tray. (2) The strength difference is caused by the inconsistency of glass fiber distribution direction and stress direction. For GFRP fiber composites, the tensile force is mainly borne by glass fiber. GFRP bars are subjected to a unidirectional tensile force, which has the same orientation as the glass fiber of GFRP bars. Therefore, the ultimate tensile strength of glass fiber can be fully played. The tray is stressed in both radial and circumferential directions. One dimensional glass fiber can not meet the stress in two directions at the same time, resulting in that only some fibers bear tensile force in each stress direction. Assuming that the glass fiber content of a tray and GFRP bar is the same, the strength of the tray under two-way stress is at least 50% lower than that of the GFRP bar under one-way stress. Accordingly, the tray strengths calculated by the distortion energy density theory and numerical simulation have reached the tray ultimate bearing capacity. However, they are about 50% of the Φ20 GFRP bar standard value of the ultimate tensile strength. The rationality that the calculated tray strength is about 50% lower than the ultimate tensile strength of GFRP bars further proves the rationality of using the theoretical calculation result of thin plate B. The distortion energy density theory considers the force and deformation when calculating the tray strength. Therefore, it is feasible and reasonable for the tray made of GFRP reinforced composites to use the distortion energy density theory to calculate the strength.

Potential Strategy to Prevent Failure of Trays
Based on the theoretical analysis and numerical simulation results of the above tray stress distribution and ultimate strength, the following improvement measures are proposed: (1) Increase the glass fiber content of the tray, and improve the strength and deformation resistance of the tray from the perspective of materials. (2) Optimize the appearance shape of the free surface of the tray and reduce the stress concentration caused by the sudden change of the thickness of the tray, so as to improve the stiffness of the tray and reduce the radial deformation of the stressed surface.

Conclusions
A stress analysis method of variable thickness tray is proposed. In this method, the variable thickness tray is simplified into an equal thickness thin plate. The stress distribution of the tray is calculated by superposition of thin plate bending theory and cavity expansion theory, and the calculation method of the tray ultimate strength is discussed.
In the nut-tray anchorage test, the failure modes are mainly nut cracking and nut sliding with GFRP bars. The tray only appears local cracking damage in the inner wall of the free surface, but the tray surface has no obvious damage, and the deformation of the tray surface is recoverable. Therefore, the elastic analysis is used in the theoretical analysis, and the linear static analysis is used in the numerical simulation. The disk load is the same as the nut-tray anchorage test, and the elastic modulus is tested by GFRP bar tensile strength test.
The thickness of the tray is the main factor affecting the size and spatial distribution of the radial deflection, radial stress, and hoop stress of the tray. In the simplified analysis, the reasonable selection of the equivalent thickness of the variable thickness tray is the key to solving the problem. By comparing the theoretical calculation and numerical simulation results, it is considered that the variable thickness tray is simplified to the equal thickness tray according to the thickness of the hole wall of the central opening of a tray, and the calculation results of tray internal force are reasonable and feasible.
The radial stress and hoop stress of the loading surface and free surface of the tray are the results of the superposition of the transverse load acting on the stressed surface of the tray and the extrusion effect of the nut on the inner wall. The stress superposition causes the tray loading surface is controlled by hoop tensile stress, and the free surface is controlled by radial compressive stress. In terms of numerical value, the hoop tensile stress on the loading surface of the tray is much greater than the radial compressive stress on the free surface.
On the loading surface of the tray, the contribution of nut extrusion to the hoop stress is dominant, accounting for 84~63%; The contribution of tray lateral load to the hoop stress occupies a secondary position, accounting for 16~37%. From the central hole wall of the tray to the edge of the tray, the hoop stress caused by the extrusion of the nut gradually decreases from large to small.
The extrusion of the nut causes the stress distribution of the central hole wall of the tray to change greatly along the thickness direction. From the loading surface of the tray to the free surface, the hoop tensile stress is dominant, and gradually transits to the radial compressive stress. The difference in the distribution of glass fiber may be the main reason for the difference in the ultimate strength between the tray and the GFRP bar. For the tray under three-dimensional stress, the distortion energy density theory must be able to theoretically evaluate the ultimate strength, which is relatively reasonable, but the tensile strength of the GFRP bar under unidirectional tension cannot be used.
The change of thickness on the free surface of a variable thickness tray will cause stress concentration. In the future, the tray design should be optimized to eliminate the influence of adverse factors.
The content of glass fiber has a great impact on the strength of the GFRP bar. GFRP bar should be tested in field application and numerical simulation to ensure that the parameter value is reasonable.
By using the distortion energy density theory to analyze the internal force of the theoretical analysis and numerical simulation calculation, it can be concluded that the tensile stress of the inner wall of the free surface of the tray is the largest, which is the weakest part of the tray, and the above analysis conclusion is the same as the damage location of the tray measured in the nut-tray anchorage test. Therefore, the theoretical analysis, numerical simulation, and nut-tray anchorage test results are mutually verified, which further proves the rationality and feasibility of the theoretical analysis method.